Visualisierung Zeckenstiche

Einführung

In dieser Aufgabe sollen die Studierenden lernen wie man Geodaten mithilfe der Packages geopandas und folium darstellen kann. Dabei werden originale Zeckenstiche und durch eine Monte Carlo simulierte Zeckenstiche dargestellt. Informationen zu den verwendeten Datensätzen finden Sie in Table 1.

Visualisierung einer Zecke. Bild von brgfx auf Freepik.

Methodik

Zuerst werden die Packages geladen und die folgende Datensätze eingelesen.

Table 1: Verwendete Datensätze

Name	Datentyp	Herkunft	Koordinatensystem	Beschreibung
zeckenstiche_original	gpkg	TickBite App	LV95	Diese Zeckenstiche wurden von 2015 bis 2020 im Umkreis von Zürich über die App TickBite erfasst. Aus dem original Datensatz wurden random 100 Zeckenstiche ausgewählt.
zeckenstiche_simuliert	gpkg	FG Geoinformatik	LV95	Simulierte Zeckenstiche aus Monte Carlo Simulation. Dabei wurden pro Zeckenstich 50 Runs mit 3 Distanzen (10m, 100m, 1000m) durchgeführt.
wald	gpkg	Swisstopo	LV95	Dieser Datensatz beinhaltet die Waldflächen und nicht-Waldflächen in und um Zürich.

Mithilfe von Folium wird eine Base Map erstellt, welche als Mittelpunkt die Stadt Zürich hat. Als es nächstes werden die Layer wald, zeckenstiche_original und zeckenstiche_simuliert hinzugeführt.

Base Map erstellen

# Create a map for all radien
## Basemap mit Koordinaten von Zürich
# Zurich latitude longitude: 47.36667, 8.55
map = folium.Map(location=[47.36667, 8.55], zoom_start = 13, tiles = "cartodbpositron")

Wald Layer hinzufügen

# Wald zur Karte hinzufügen
is_wald = wald[wald["Wald_text"] == "ja"]
is_wald.to_crs(4326).explore(
    m = map,           # pass the map object
    color = "green",   # use green color on all points
    tooltip = False,   # hide tooltip
    highlight = False, # disable highlighting
    name = "Wald"      # name of the layer in the map
)

Original Zeckenstiche hinzufügen

# Orginal Zeckenstiche zur Karte hinzufügen
zeckenstiche_original_gpd.to_crs(4326).explore(
    m = map,                                     # pass the map object
    color = "red",                               # use red color on all points
    marker_kwds = dict(radius = 1, fill = True), # make marker radius 1px with fill
    tooltip = "ID",                              # show "ID" column in the tooltip
    name = "Original Zeckenstiche"               # name of the layer in the map
)

Simulierte Zeckenstiche hinzufügen

# Visualisieren der Simulierten Zeckenstiche
zeckenstiche_sim_gpd.to_crs(4326).explore(
    m = map,
    column = "Radius",
    cmap = "viridis",
    legend = True,
    legend_kwds = dict(colorbar = False),
    marker_kwds = dict(radius = 1, fill = True),
    tooltip = ["ID", "Radius", "Run"],
    name = "Simulation Zeckenstiche",
)

Layer Control hinzufügen

folium.LayerControl().add_to(map)

map

Make this Notebook Trusted to load map: File -> Trust Notebook